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DESCRIPTION 

TITLE OF THE INVENTION 

Methods of Designing Optimal PID Controllers 

TECHNICAL FIELD OF THE INVENTION 

This invention relates to the design of the structure of a multivariable PID controller 
and the optima! choice of its PID parameters. 

BACKGROUND OF THE INVENTION 

A traditional PID controller is used to control an industrial process. The process 
variable (PV) goes into the PID controller, which calculates the controller output (CO) 
according to a PID control equation. This CO is then converted to an analog signal, 
which is sent to the process so that the said PV can track a user specified value called 
set point (SP). The said SP can change with time. The performance of a PID controller 
depends on the choice of its three PID parameters. For independent form of PID 
controllers these three PID parameters are the proportional gain Kp, the integral gain 
Ki, and the derivative gain Kd. For dependent form of PID controllers these three PID 
parameters are the gain K, integral time Ti, and derivative time Td. In traditional PID 
controllers the said PV, SP, CO, and PID parameters are all scalars. We call this kind 
of PID controllers the single-input single-output (SISO) PID controllers. The Ziegler- 
Nichols PID controller tuning method is the major one of the many methods for 
finding the values of PID parameters. 

DETAILED DESCRIPTION OF THE INVENTION 

In this invention the SISO PID controller is extended to the multiple-input multiple- 
output (M1MO) PID controller that has n process variables PV1, PV2, and PVn 
and m controller outputs CO 1 , C02, . . . , and COm, where m and n are positive 
integers. Corresponding to PV1, PV2, and PVn there are n set points SP1, SP2, 
and SPn. In this case PV becomes a vector with PV1, PV2, .. ., and PVn being its first, 
second, and n-th component, CO becomes a vector with COl, C02, and COm 
being its first, second, . .., and m-th component, SP becomes a vector with SP1, 
SP2, and SPn being its first, second, . . ., and n-th component, and the PID 
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control equation becomes CO(k) = CO(k-l) + Kl*SP(k)*T + Kl*a(k,l) + K2*a(k,2) 
+ . . . + Kj*a(k,j), where k is the discrete time, T is the sampling period, j is a positive 
integer, Kl, K2, Kj are m by n PED parameters, a(k,l) = [-PV(k)]*T, and a(kj) = 
[a(k,j-l) - a(k-l ,j-l)]/T for j > or = 2. It is important to note that 

1 . an MIMO PID controller is able to take into account the interaction among the n 
process variables and m controller outputs, which can not be achieved by simply 
applying S1SO PID controllers to each of the n process variables, and 

2. there is no set point in any of a(k, 1), a(k,2), . . . , and a(kj), which can avoid the 
unwanted sudden change in CO when SP changes with time. 

The next problem of designing the optimal PID controller is to find the best values for 
the PID parameters Kl, K2, and Kj. An optimization based method for solving 
this problem consists of the following four steps: 

1 . Convert the PID control equation into discrete time form if it is not in discrete 
time form. 

2. Build a discrete time linear model for the process that is to be controlled by the 
said PID controller. 

3. Form the discrete time closed loop transfer function from said vector SP to said 
vector PV. 

4. Find the best PID parameters by using an optimization algorithm which minimizes 
the largest modulus of all poles of the discrete time closed loop transfer function 
obtained at step 3, where the modulus of a pole is defined to be the absolute value 
of the complex number which represents the pole. If the PID parameters are 
subject to some constraints, then a constrained optimization algorithm can be used 
which minimizes the largest modulus of all poles of the discrete time closed loop 
transfer function obtained at step 3 and at the same time guarantees that all user 
prescribed constraints on the PID parameters are satisfied. 

PID controllers with their parameters so obtained guarantee that PV can track SP 
quickly. 



